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In this paper, the optimal incorporation of distributed generation (OIDG) 
units for the maximization of the system’s loadability is investigated (based 
on sizing and siting). To this end a new computational approach for 
computing maximum loadability of the system is developed. This approach 
has been compared with the classical one on different radial test systems 
(RTS) and is found to be faster and more accurate. The OIDG problem is 
formulated mathematically as an optimization problem with the objective 
function to maximize system’s loadability, the imposed constraints are; 
voltage limits, thermal limits and DG penetration level. The optimization 
algorithm used is to solve the OIDG problem is the Vortex searching 
algorithm (VS). The tested radial distribution systems are the standard 33-bus 
and 69-bus systems. This paper also discusses some other interesting findings 
about VS algorithm. 
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1. INTRODUCTION 

The word distributed generation (DG) is subjective in nature and its usage changes from region to 

another. For instance, in Anglo-American countries the term ‘embedded generation’ is used, however North 
American countries prefer it to call as ‘dispersed generation’, whereas Europe and parts of Asia use the 
terminology of ‘decentralized generation’ [1]. 

The definition of DG depends on the purpose like the size (MW), voltage level (kV), technology, 
cogeneration and others [1] defines. At present, council on large electric systems (C1GRE) considers DG to 
be not centrally planned, nor centrally dispatched and usually connected to the distribution network, smaller 
than 50-100 MW [2]. Various authors describe DG as a power generation source incorporated in the 
distribution network ranging from a few kW to a few tens of MW [1]. Some example of DG categories and 
ratings are micro-distributed generation (~1 W<5 kW), small distributed generation (5 kW<5 MW), medium 
distributed generation (5 MW<50 MW) and large distributed generation 50 MW<300 MW. 

A DG can be a traditional generator - based on combustion of diesel reciprocating generator and 
natural gas-turbine-and/or a non-traditional generator - relying on storage devices, fuel cells, or renewable 
energy source (such as wind turbine and photovoltaic) [1-3]. Advantages of DG over centralized generation 
are numerous including [4-7]: 
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a. Power loss reduction 

b. Voltage profile improvement 

c. Enhanced system stability 

d. Reduced pollutant emission 

e. Relieved transmission and distribution lines. 

However, there are several issues related with DG incorporation (i.e. placement and sizing) in the 
existing networks. The incorporation of DG alters the nature of network from passive one to that of active 
one, which has a significant impact in network operations. Non-optimized incorporation of DG can result in 
deterioration of the network performance. Therefore, the optimal incorporation of distributed generation 
(OIDG) is of paramount importance. 

The problem of OIDG has been solved using analytical approaches. The initial work in DG 
placement was carried out by [8] based on the famous 2/3 rule of thumb. This work has been extended in [9] 
for three different kinds of loads. Later on, other analytical approaches have also been proposed to determine 
optimum DG placement as well as DG size. In [2], optimum size of DG is determined by using loss 
sensitivity based analytical equation, whereas, optimum DG location is determined by exact loss equation. 
Reference [10] covers the equivalent current injection based calculation of loss sensitivity factor injection by 
means of two matrixes: Bus-Injection to Branch-Current (BIBC) matrix and Branch-Current to Bus-Voltage 
(BCBV) matrix. A conventional iterative search algorithm is presented in [11] which deals with OIDG based 
on loss reduction and cost of DG size. Although this method is simple, it consumes a lot of time. Reference 
[12] ntroduces an improvement in analytical method for identification of optimal location for DG installation 
along with optimal power factor to achieve loss reduction in distribution networks. 

Algorithms, based on multi-objectives optimization, have also been proposed for optimum DG 
placement and sizing. Methods based on Genetic Algorithm (GA) are proposed in [13-14] for optimizing DG 
placement and sizing, with various objective functions. Reference [15] has combined GA with Particle 
Swarm Optimization fPSO) for OIDG, solving multi-objective functions that include power losses, voltage 
regulation and voltage stability. GA-Fuzzy approach is used in [16] to optimize location of DG. Objective 
function also covers reduction in power system losses and system loading as well as the profit maximization 
for Distribution Companies. A dynamic programming approach is introduced in [17] to ascertain an 
optimized DG location that would ensure maximum profit. Reference [18] has used Ant Colony Search 
(ACS) Algorithm for reliability based OIDG location wise. In [19] PSO algorithm is implemented for 
optimum placement of DG while minimizing electricity cost for consumers. Reference [20] has used cuckoo 
search algorithm for optimum DG placement to obtain better voltage profile and also to decrease the power 
losses of the system. 

Various other researchers have also considered OIDG as a multi-objective function optimization 
problem. Their fitness functions may include improved voltage regulation, reduction in power losses, and 
voltage stability [15]. Reference [13] has considered cost benefit due to reduction of power losses in case of 
DG placement in fitness function and ensure that voltage profile and reliability of the system is improved. It 
has been reported in [21] that incorporation of DGs gives rise to the loadability of overall 
distribution network. 

In this paper, the optimal simultaneous incorporation of multiple DGs based on maximization of 
system loadability is investigated where the maximum loadability is determined using a new fast and 
accurate approach. The optimization algorithm used to solve the OIDG is the Vortex Searching algorithm 
(VS). This remainder of this paper is structured as follows. Section 2 presents the formulation of the OIDG 
problem. The optimization algorithm used in this paper, VS, is explained in section 3. Section 4 includes the 
discussion of main results obtained in this paper. Findings of this paper are concluded in section 5. 


2. FORMULATION OF THE OIDG PROBLEM 

Explaining The primary purpose of OIDG problem is to determine the optimal placement and size 
of DG for improved system performance. The OIDG problem can be formulated as follows [22]: 


Minimize 

/(x) 

(1) 

Subject to 

< 7 j(x) = 0 i = 1: k 

(2) 

and 

hj(x) <0 j = l:m 

(3) 

where 

X = [x^x- 2 ,... ,x n ] 

(4) 
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where: 

/(x) is the objective function. 

x is the vector of design variables, x- L is the ith design variable and n represents the number of 
design variables. 

g(x) denotes the set of equality constraints and k is the number of equality constraints. 
h(x) is the set of inequality constraints and m is the number of inequality constraints. 

Design variables 

As previously mentioned, this paper aims to ascertain the optimal size and location of DGs to be 
inserted in the system. Therefore, the following design variables are considered. 

Size 

The size of DG means the amount of power generated by the DG unit. In following sections, the 
active power of DG, denoted by P DG , is considered as a design variable and the reactive power of DG, Q DG , is 
calculated using the following expression: 

Qdg = Pdg x tan(cos _1 (PF)) (5) 

The size of DG has to respect the following constraint: 

0 — 2fc=i P D G k — P L oad (6) 

where: PF stands for power factor. 

Location 

Location of DG means the optimal bus to insert the DG unit. The location or position of DG has to 
respect the following constraint: 

2 < DG location < nb (7) 

where: nb is the number of buses. 

Furthermore, if there are more than one of DG unit, the following constraint has to be 
respected also: 

DGlocation 1 T DG[ ocat j on 2 ^ •" ^ DGi ocat j on n dg (8) 

Power factor 

In some of the investigated cases of this paper, the PF of the DG unit is imposed and in some of 
them the PF is considered as a design variable and it has to be determined. 

Objective function 

Here, the main objective function is system maximum loadability. 

Maximum loadability 

Different definitions and different terminologies of loadability and maximum loadability have been 
reported in [23-28]. These definitions vary in terms of constraints, stability of system and the representation 
of loadability. 

Maximum loadability or maximum loading margin is used to refer maximum power load that could 
load the system to the point of instability from the base load. Maximum loadability of the system indicates 
how much the system must be stressed in order to reach the instability. Voltage instability is considered a 
crucial parameter throughout this paper, thus the maximum loadability of the system can be defined as the 
maximum increase in power system load from the base load till the point where voltage collapses. Maximum 
loadability of the system represents the loading factor in multiple of base case load till the point of voltage 
collapse and represented with the variable A, max . The maximum loadability of the system can also be referred 
as or Voltage Stability Margin (VSM) or Voltage Stability Limit (VSL). 

Proposed approach 

In this paper, a new approach for calculating A max is proposed. This approach is much faster than 
the classical one proposed and implemented in [21], The flowchart of the proposed approach is given in 
Figure 1. 
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Figure 1. Flowchart for the calculation of A max . 


To find the maximum loadability i.e. A max , on each bus, both active and reactive power are 
increased using the following: 

P — P 0 X A (9) 

- Qo X A (10) 
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where: P 0 and Q 0 are the initial active and reactive powers of combined loads, respectively. Similarly, P and 
Q represent the new active and reactive powers (loads), respectively and A is the loading factor. 

However, instead of keep increasing A using a smal step as in the classical approach [21], here we 
start from an initial interval [a, b], then A is calculatd as the midpoint of this interval. After that, P and Q are 
increased using equations (9) and (10). The power flow equations are solved, if it is successful the new 
interval will be [a. A] otherwise the new interval will be [A, b]. This process is continued until the interval [a, 
b] becomes smaller than a predefined tolerance (£). 

Comparative study 

In order to show the speed and efficiency of the proposed approach compared to the one proposed in 
[21], several tests have been achieved as shown in Table 1. In this table, the A max and the time taken to 
compute the A max are reported. 

It can be noticed from Table 1 that, the A max is the same using both approaches for all cases, 
however, there is a huge time difference between both approaches. For instance for the 10-bus radial system 
the time has been reduced by 94.64%, for the 69-bus radial system the time has been reduced by 99.14% and 
for the 118-bus radial system the time has been reduced by 99.07% using the proposed approach compared to 
the one proposed in [21]. 

Moreover, it is worth mentioning that, this time has been computed for one simulation, therefore, if 
this simulation is repeated a thousand times as when solving the OIDG problem the gain in time becomes 
extremely important. 


Table 1. Comparison of the Proposed Approach with the Classical One 


Test system 

Classical approach 

Proposed approach 

^max 

Time (s) 

^max 

Time (s) 

10-bus radial system [29] 

2.064 

1.698 

2.063 

0.091 

16-bus radial system [21] 

7.550 

7.852 

7.549 

0.110 

33-bus radial system [30] 

3.407 

9.648 

3.406 

0.125 

69-bus radial system [31] 

3.210 

28.813 

3.210 

0.245 

85-bus radial system [32] 

2.599 

31.958 

2.599 

0.326 

118-bus radial system [33] 

2.465 

46.933 

2.464 

0.438 


Constraints 

For OIDG, both equality and inequality constraints are present in problem formulation. 

Equality constraints 

Pgi =Pi+V t VjYtj cos (St - Sj - fly) (11) 

Qgt = Qi + Vi VjYij sin (Aj - Sj - fl y ) (12) 

where; 

rij = 1; 2; .n n 

P g iis generator active power output at bustij 

Qgiis generator reactive power output at bustij 

Pi is active power demand at bus ni 

Qi is reactive power demand at bus rij 

Vi is voltage of bus ni ; S t is phase voltage angle at bus ni; 

N — ni _ 1 is total number of branches in the given radial distribution system. 

N represents total number of buses in given radial distribution system 

Inequality constraints 
Voltage limits 

The nodal voltages must be constrained between their minimum and maximum limits. Therefore, 
voltage limits constraint can be formulated as follows: 

F min < ^ < V max (13) 

Thermal limit 

Power flow through any distribution feeder must comply with the thermal capacity of the line. 
Therefore, the thermal limit constraints can be formulated as follows: 
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|S t l<|S i max U = l. N 


(14) 


where: 5) is the apparent power at branch i and S ™ ax is the maximum apparent power at branch i. 


Maximum capacity 

The maximum capacity of DGs that can be installed in the network is restricted by DG penetration 
level ( DG pl ). Therefore, there is a capacity constraint at unity power factor of DG generating power that is 
given by following expression: 


(DG PL XP load ) 

n DG 


(15) 


If both active and reactive powers are being generated by DG, then the precedent equation can be 
modified as: 




(DG pl xSiq ad ) 
n DG 


(16) 


where N DG is the total number of DG units incorporated in the power system; DG PL represents maximum 
penetration level as a percentage of peak load of the system; P ioad and S load are denoting total active power 
load and total apparent power load of the system, respectively. 


Constraints handling procedure 

It is worth to mention that the inequality constraints are taken into account using the penalty 
method. In this method a penalty term is added to the primary objective function which consists of 
multiplying a penalty parameter by a penalty factor. The penalty parameter is a measure of how much the 
constraints are violated, if there is no violation this parameter is null. This method has been implemented and 
used by the VS algorithm to handle inequality constraints. 

Evaluation indices 

The evaluation of the OIDG on system performance is performed using the following indices where 
the subscript (0) represents values before the incorporation of DG and the subscript DG represents the values 
calculated after the incorporation of DG [21], 


Active losses reduction index 

The Active Losses Reduction f ALR) index is be given by: 


— P Losses 0 PLosses DG n/ 
^LossesQ 


(17) 


where: The active power losses are set as follows: 

bosses =52=ifl k |/* 2 1 (18) 

where: R k is the resistance of branch k, I k represents the current passing through branch k and N denotes the 
total number of branches. 

Reactive losses reduction index 

The Reactive Losses Reduction (RLR) index is given by: 


_ Qrossesp QLossespg (19) 

QLossesp 

System loadability improvement index 

The system loadability improvement (SLI) index is given by: 

SL] = W DG -W 0 % (20) 

^maxo 
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Voltage profile improvement index 

The voltage profile improvement (VPI) index is given by: 

VPI = vp °— VPp G o/ 0 (21) 

VPo 

Where the voltage profile is defined as follows: 

VP = - V rated ) 2 (22) 


Where Prated is the rated voltage (taken as 1 p.u. in this paper). 

Voltage stability enhancement index 

The voltage stability enhancement index (VSIM) is as follows: 

VSIM = VSlpG ~ VS1 ° % (23) 

vsi 0 ' 

and voltage stability index (VSI) will be expressed as [34]: 

VSI = min n t = 2,3,..., nb (24) 

where SI of node rij is given by: 

SI(n,) = \V mi \ 4 - 4 [PniindRm + Qm(ndX ni ]\V ml \ 4 - 4[P ni {n{)X ni + Q ni (ndR ni \ 2 (25) 

where V mi is the voltage of bus mi, P nl (nj) is total real power load fed through bus ni, Q n i{nd is 
total reactive power load fed through bus ni, R ni is the resistance of branch I, X ni is the reactance of branch i. 

Number of buses violating voltage limit 

Number of buses violating voltage limit (NBVV), as indicated by its name, counts the numbers of 
busses where voltage limits (upper of lower limits) are not respected. 


3. OPTIMIZATION ALGORITHM 

The VS algorithm is a single-solution based metaheuristic that evolves to optimality using a 
generation and a replacement procedure. In the generation phase N candidate solutions are generated inside a 
radius around the best solution using a specified procedure and, in the replacement, phase a solution replaces 
the current best solution. This process is repeated until a termination criterion is met. It is reported in [35] 
that the VS has excellent performance for the optimization of mathematical functions. In this paper is used as 
the optimization algorithm for solving the OIDG problem. 

Algorithm description 

The flowchart of the VS algorithm is sketched in Figure 1. The main steps of the VS algorithm are 
explained below: 

Step 1: in the initialisation phase, first the centre of the search space is selected as the first generated 
solution. Then, the initial radius is computed using the following expression: 

r 0 = a 0 -gammaincinv(jc,a 0 ) (26) 

where: gammaincinv is the inverse incomplete gamma function and <r 0 is the initial standard 
distribution calculated as follows: 

max(t/Z?)-min(ZJ?) 

Oo =--- (17) 

where: UB is the upper bound of design variables, LB is the lower bound of the design variables, x 
is a random variable and a 0 is the shape parameter selected as 1. 
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After that, the best solution and fitness obtained so far referred to as gb es t and S bes t, respectively are 
initialized with the first solution generated. 

Step 2: in this step N candidate solutions (with N the size of solutions) are generated in the 
neighborhood of g bes t using a Gaussian distribution where this neighborhood is characterized with the 
following radius: 


r = a 0 -gammaincinv(x, a t ) 


(28) 


where the shape parameter a t is decreased using the following expression: 


&t — a q 


t 

Maxltr 


(29) 


where t is iteration number and Maxltr is the maximum number of iterations. 

It is worth to mention here that; this radius is of paramount importance in the convergence of the VS 
algorithm. For a given optimization problem, this radius is initialized with a high value (as in equation (1)) in 
order to promote the exploration phase in the first iterations and then it is decreased along with the number of 
iterations increases in order to promote the exploitation phase. It is reported in [35] that until a certain 
number of iterations is reached, the radius varies approximately linearly and then it decreases significantly. 
This behavior allows to have a good balance between exploration and exploitation phases. 

Step 3: in the previous step some of the generated solutions may lay outside the search space, 
therefore they have to be brought back inside the search space using the following expression: 

(LB 1 + rand ■ (UB l - LB 1 ) if S l k < LB‘ 

Si = ] S l k if LB 1 < Si < UB l (30) 

[LB 1 + rand ■ (UB l - LB 1 ) if 4 > UB l 

Step 4: in this step the generated solutions are evaluated in order to calculate their fitness’s, the best 
solution among them is identified and noted as gbest_New and Sb es t_New- Then, if Sbest.New is found to be better 
than S beS b therefore g bes , and S bes t are updated. 

Step 5: in this step the radius is updated using equation (28). 

Step 6: finally, if the number of iterations exceeds Maxltr the process is stopped and g bes t and S best 
are displayed otherwise the process is repeated from Step 2. 

The VS algorithm as it is proposed in is [35] does not handle discrete or integer design variables. 
Therefore, for the purpose of the OIDG problem we have added some modifications on the initial version of 
the VS algorithm to cope with this issue. 


4. RESULT AND DISCUSSION 

Provide In this paper, 62 case studies have been investigated using 2 test systems. More details 
about test systems and investigated cases are given in the following two subsections. 

Test systems 

In this paper, two test systems are considered: 33-bus and 69-bus RTS. The main characteristics of 
these systems are given in Table 2. Moreover, the single line diagrams of first and second test systems are 
shown in Figure 3 and Figure 4, respectively. 


Table 2. The Main Characteristics of the Investigated Test Systems 


System 

3 3-bus RTS 

69-bus RTS 

System Characteristics 

Data is given in [30] 

Data is given in [31] 

Buses 

33 

69 

Branches 

32 

68 

Total active load (MW) 

3.7150 

3.8021 

Total reactive load (Mvar) 

2.3000 

2.6945 

Ploss (MW) 

0.2110 

0.2250 

Qloss (Mvar) 

0.1430 

0.1021 

^■max 

3.4065 

3.2102 

VSI 

1 

1 

VP 

0.1338 

0.0993 

■ymax 

1.1 

1.1 

ymin 

0.95 

0.95 

S m “ (MV A) 

5 

5 
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Figure 2. Flowchart of the vs algorithm 
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23 24 2S 



Figure 3. Single line diagram of the 33-bus RTS 


J6 J? J4 4J <42 41 42 4J 44 O * 



DG types 

Mainly, there are four different categories of DGs considered in literature as reported in [24]. 
Category 1: DG supplies active power only, e.g. solar system. 

Category 2: DG supplies reactive power only, e.g. synchronous compensators. 

Category 3: DG supplies active power but consumes reactive power, e.g. induction generator. 
Category 4: DG supplies both active and reactive power e.g. synchronous generator. 

In this work, where the PF is considered as unity the DG units are considered from Type 1 whilst 
when the PF is considered as a design variable the DG units are considered as active-reactive power sources 
(i.e. Type 4). 

Investigated cases 

As previously mentioned, in this paper, 62 case studies are investigated as shown in Table 3. 
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CASE 1.1 through CASE 1.31 are related to the first test system whilst CASE 2.1 trough CASE 
2.31 are related to the second test system. CASE 1.1 and CASE 2.1 are base cases where no DG is inserted, 
these cases are investigated to serve as a reference for the remaining cases for comparison purposes. For the 
remaining cases the objective is to improve A max . 

CASE 1.2 to CASE 1.6 are related to one DG unit using a unity power factor and different 
penetration levels. CASE 1.7 to CASE 1.11 are related to one DG unit with different penetration levels but 
with an optimized power factor. 

CASE 1.12 to CASE 1.16 are related to two DG units using a unity power factor and different 
penetration levels. CASE 1.17 to CASE 1.21 are related to two DG units with different penetration levels but 
with an optimized power factor. 

CASE 1.22 to CASE 1.26 are related to three DG units using a unity power factor and different 
penetration levels. CASE 1.27 to CASE 1.31 are related to three DG units with different penetration levels 
but with an optimized power factor. 

Finally, the same description can be made for the cases that are related to the second test system. 


Table 3. Summary of the Studied Cases 


CASE 

No. DG 

DGpl 

PF 

Test system (radial system) 

CASE1.1/CASE2.1 

- 

- 

- 

33-bus/69-bus 

CASE1.2/CASE2.2 

1 

0.25 

1 

33-bus/69-bus 

CASE1.3/CASE2.3 

1 

0.5 

1 

33-bus/69-bus 

CASE1.4/CASE2.4 

1 

0.75 

1 

33-bus/69-bus 

CASE 1.5/CASE 2.5 

1 

1 

1 

33-bus/69-bus 

CASE 1.6/CASE 2.6 

1 

2 

1 

33-bus/69-bus 

CASE 1.7/CASE 2.7 

1 

0.25 

Design variable 

33-bus/69-bus 

CASE1.8/CASE 2.8 

1 

0.5 

Design variable 

33-bus/69-bus 

CASE 1.9/CASE 2.9 

1 

0.75 

Design variable 

33-bus/69-bus 

CASE1.10/CASE 2.10 

1 

1 

Design variable 

33-bus/69-bus 

CASE1.11/CASE 2.11 

1 

2 

Design variable 

33-bus/69-bus 

CASE1.12/CASE 2.12 

2 

0.25 

1 

33-bus/69-bus 

CASE1.13/CASE 2.13 

2 

0.5 

1 

33-bus/69-bus 

CASE1.14/CASE 2.14 

2 

0.75 

1 

33-bus/69-bus 

CASE1.15/CASE 2.15 

2 

1 

1 

33-bus/69-bus 

CASE1.16/CASE 2.16 

2 

2 

1 

33-bus/69-bus 

CASE1.17/CASE 2.17 

2 

0.25 

Design variable 

33-bus/69-bus 

CASE 1.18/CASE 2.18 

2 

0.5 

Design variable 

33-bus/69-bus 

CASE1.19/CASE 2.19 

2 

0.75 

Design variable 

33-bus/69-bus 

CASE1.20/CASE 2.20 

2 

1 

Design variable 

33-bus/69-bus 

CASE1.21/CASE 2.21 

2 

2 

Design variable 

33-bus/69-bus 

CASE1.22/CASE 2.22 

3 

0.25 

1 

33-bus/69-bus 

CASE1.23/CASE 2.23 

3 

0.5 

1 

33-bus/69-bus 

CASE1.24/CASE 2.24 

3 

0.75 

1 

33-bus/69-bus 

CASE1.25/CASE 2.25 

3 

1 

1 

33-bus/69-bus 

CASE1.26/CASE 2.26 

3 

2 

1 

33-bus/69-bus 

CASE1.27/CASE 2.27 

3 

0.25 

Design variable 

33-bus/69-bus 

CASE1.28/CASE 2.28 

3 

0.5 

Design variable 

33-bus/69-bus 

CASE1.29/CASE 2.29 

3 

0.75 

Design variable 

33-bus /69-bus 

CASE1.30/CASE 2.30 

3 

1 

Design variable 

33-bus/69-bus 

CASE1.31/CASE 2.31 

3 

2 

Design variable 

33-bus/69-bus 


The developed program has been implemented using the commercial MATLAB software (version 
R2015a). The simulation runs were performed using the proposed VS based approach with n = 50, and a 
maximum of 150 iterations. 

Results for the 33-bus radial test system 
Results for one DG unit incorporation 

The obtained results for the incorporation of one DG unit for the 33-bus RTS are presented in 
Table 4. The following comments can be made on this table: 

a. For the base case without DG (i.e. CASE 1.1), T max = 3.406. In this case NBVV=21. 

b. For CASE 1.2, CASE 1.3 and CASE 1.7, NBVV is equal to 6, 4 and 5 respectively. Even though the 
NBVV has been reduced from 21 in the base case to 6, 4 and 5 in these cases, the network is still not 
fully optimized due to incorporation of only one DG and due to the limitation on the penetration level. 

c. The location of the DG unit is between buses 8 and 13. 
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d. For CASE 1.4, using only one DG unit but with a penetration level of 75% the network is optimized and 
the NBVV=0. The maximum loadability in this case has been increased from 3.406 to 4.195 compared 
with the base case i.e. SLI=23.14%. 

e. For CASE 1.5, the SLI is 25.07% while it is 26.88% in CASE 1.6 where no constraint is imposed on the 
penetration level of DG. 

f. For CASE 1.8, even though only one DG unit is incorporated using a penetration lever of 50% the 
network is optimized (NBVV=0) compared with CASE 1.3 because the PF has been computed as 0.88 
instead of 1 for CASE 1.3. 

g. For CASE 1.9 and 1.10 the SLIs are 25.93% and 26.88% compared with the base case. 

h. The PV curves of CASE 1.1, CASE 1.5 and CASE 1.10 are displayed in Figure 5. 

i. The voltage profile is improved in all cases in comparison to the base case as shown by the VPI. The 
best improvement is obtained for CASE 1.6 where VP1=88.206%. The voltage profile of CASE 1.1, 
CASE 1.5 and CASE 1.10 are displayed in Figure 6. 

j. The active power losses are reduced in some cases as for CASE 1.8 and CASE 1.9 and they are 
increased in other cases as for CASE 1.4, CASE 1.5, CASE 1.6, CASE 1.10 and CASE 1.11 

k. The evolution of the objective function for CASE 1.5 and CASE 1.10 is given in Figure 7. 


Results for two DG units incorporation 

The obtained results for the incorporation of two DG units for the 33-bus RTS are given in Table 5. 

The following comments can be made on this table: 

a. When incorporating two DG units only one case is not optimized where NBVV=15 that is CASE 1.12. 

b. The location of the first DG unit is between buses 14 and 18 while the location of the second DG unit is 
between buses 31 and 32 for all cases. 

c. The most important SLI is for CASE 1.21 where it is equal to 47.86%. An illustration of the PV curves 
of CASE 1.1, CASE 1.15 and CASE 1.20 is given in Figure 8. 

d. The voltage profile is improved in all cases compared with the base case. An illustration of VP for 
CASE 1.1, CASE 1.15 and CASE 1.20 is given in Figure 9. 

e. Active power losses are reduced in almost all cases except CASE 1.16 where there is an increase of 
48.06%. 

f. When the PF is not fixed and it is considered as a design variable, it varies approximatively between 
0.82 and 0.95. 

g. The evolution of the objective function for CASE 1.15 and CASE 1.20 is sketched in Figure 10. 

Results for three DG units incorporation 

The obtained results for the incorporation of three DG units for the 33-bus RTS are given in Table 6. 

The following comments can be made on this table: 

a. When incorporating three DG units only one case is not optimized where NBVV=5 that is CASE 1.22. 

b. The location of DG units is mainly at busses 10, 14, 15, 16, 17, 18, 30, 31 and 32. 

c. The best SLI obtained is for CASE 1.31 where it is equal to 48.64%. An illustration of the PV curves of 
CASE 1.1, CASE 1.25 and CASE 1.30 is given in Figure 11. 

d. The voltage profile is improved in all cases compared with the base case. An illustration of VP for 
CASE 1.1, CASE 1.25 and CASE 1.30 is given in Figure 12. 

e. Active power losses are reduced in all cases compared with the base case except for CASE 1.26. 

f. When the PF is not fixed and it is considered as a design variable, it varies approximatively between 
0.81 and 0.90. 

The evolution of the objective function for CASE 1.25 and CASE 1.30 is given in Figure 13. 
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Table 4. Obtained Results for the 33-bus RTS for One DG Unit 


Location Size (P [MW]) Performance evaluation indices 


CASE# 

DG, 

dg 2 

DG 

DG, 

dg 2 

DG 3 

PF 

DGpl 

4max 

P Losses 

VSI 

VP 

ALR 

RLR 

SLI 

VPI 

VSIM 

NBVV 

CASE 

1.1 

- 

- 

- 

- 

- 

- 

- 

- 

3.40 

6 

0.211 

1 

0.134 

- 

- 

- 

- 

- 

21 

CASE 

1.2 

13 

- 

- 

0.930 

- 

- 

1 

25 

3.80 

5 

0.132 

1 

0.052 

37.23 

1 

38.0 

81 

11.6 

96 

61.3 

71 

0 

6 

CASE 

1.3 

8 

- 

- 

1.858 

- 

- 

1 

50 

3.81 

4 

0.118 

1 

0.034 

43.95 

7 

41.6 

60 

11.9 

54 

74.5 

55 

0 

4 

CASE 

1.4 

13 

- 

- 

2.780 

- 

- 

1 

75 

4.19 

5 

0.248 

0.82 

7 

0.027 

17.69 

3 

26.3 

55 

23.1 

35 

79.8 

05 

17.275 

0 

CASE 

1.5 

10 

- 

- 

3.499 

- 

- 

1 

100 

4.26 

0 

0.275 

0.82 

5 

0.025 

30.48 

3 

47.7 

43 

25.0 

70 

81.6 

18 

17.517 

0 

CASE 

1.6 

8 

- 

- 

4.809 

- 

- 

1 

- 

4.32 

2 

0.327 

0.83 

3 

0.016 

54.80 

9 

87.9 

17 

26.8 

76 

88.2 

06 

16.678 

0 

CASE 

1.7 

12 

- 

- 

0.922 

- 

- 

0.84 

4 

25 

3.86 

2 

0.105 

1 

0.036 

50.18 

9 

50.8 

85 

13.3 

73 

73.1 

20 

0 

5 

CASE 

1.8 

13 

- 

- 

1.932 

- 

- 

0.88 

4 

50 

4.19 

8 

0.131 

0.82 

6 

0.028 

38.14 

5 

34.2 

53 

23.2 

21 

79.3 

27 

17.416 

0 

CASE 

1.9 

9 

- 

- 

3.099 

- 

- 

0.94 

7 

75 

4.29 

0 

0.155 

0.82 

4 

0.021 

26.69 

3 

13.4 

04 

25.9 

30 

84.4 

84 

17.571 

0 

CASE 

1.10 

8 

- 

- 

4.200 

- 

- 

0.98 

8 

100 

4.32 

2 

0.222 

0.83 

1 

0.016 

5.393 

29.7 

02 

26.8 

76 

88.1 

54 

16.944 

0 

CASE 

1.11 

8 

- 

- 

4.522 

- 

- 

0.99 

8 

- 

4.32 

3 

0.273 

0.83 

2 

0.016 

29.27 

1 

57.8 

50 

26.9 

19 

88.1 

38 

16.845 

0 





Table 5. 

Obtained Results for the 33-bus RTS for Two DG Units 







CASE 


Location 


Size (P [MW]) 








Performance evaluation indices 


# 

DGi 

dg 2 

dg 3 

DGi 

DG 2 DG 3 

PF 

DGpl 

^max 

P Losses 

VSI 

VP 

ALR 

RLR 

SLI 

VPI 

VSIM 

NBVV 

CASE 








3.406 

0.211 

1 

0.1 






21 

- 

- 

- 

- 

- 

- 

_ 

- 

- 

- 

- 

. 

1.1 










34 







CASE 

18 

32 


0.46 

0.46 

1 

25 

3.828 

0.119 

1 

0.0 

43.4 

44.1 

12.3 

58.2 

0 

15 

1.12 


6 

6 

56 

16 

01 

84 

14 

CASE 

17 

32 


0.92 

0.92 

1 

50 

4.176 

0.100 

1 

0.0 

52.7 

49.5 

22.5 

87.8 

0 

0 

_ 


1.13 


7 

8 

16 

53 

81 

76 

88 

CASE 

17 

32 


1.39 

1.39 

1 

75 

4.474 

0.140 

0.91 

0.0 

33.7 

22.8 

31.3 

96.9 

8.732 

0 

_ 


1.14 


3 

2 

3 

04 

67 

96 

48 

32 

CASE 

16 

32 


1.85 

1.85 

1 

100 

4.748 

0.209 

0.83 

0.0 

0.94 

11.3 

39.3 

90.7 

16.687 

0 

1.15 


7 

7 


3 

12 

0 

24 

89 

61 

CASE 

14 

31 


1.89 

2.93 

1 


4.959 

0.312 

0.82 

0.0 

- 

- 

45.5 

78.4 

17.695 

0 

_ 


_ 

48.0 

64.8 

1.16 



7 

4 



3 

29 

59 

42 

82 

23 

CASE 

17 

32 


0.46 

0.46 

0.8 

25 

3.991 

0.082 

1 

0.0 

61.3 

62.1 

17.1 

75.0 

0 

0 

1.17 

" 

4 

5 

52 

33 

63 

03 

58 

58 

CASE 

17 

32 


0.89 

0.89 

0.8 

50 

4.505 

0.051 

0.89 

0.0 

75.9 

71.3 

32.2 

96.3 

10.148 

0 

1.18 


9 

9 

24 

9 

05 

42 

38 

51 

85 

CASE 

14 

32 


1.43 

1.43 

0.8 

75 

4.920 

0.071 

0.82 

0.0 

66.4 

59.8 

44.4 

85.4 

17.683 

0 

1.19 


3 

3 

75 

3 

19 

22 

26 

21 

46 

CASE 

14 

32 


1.38 

1.97 

0.9 

100 

5.024 

0.099 

0.82 

0.0 

53.3 

43.4 

47.4 

79.8 

17.713 

0 

_ 


74 

1.20 


5 

1 

02 

3 

27 

02 

12 

37 

CASE 

15 

31 


1.38 

2.30 

0.9 


5.037 

0.124 

0.82 

0.0 

41.1 

31.2 

47.8 

79.5 

17.706 

0 

1.21 

" 

1 

5 

45 

" 

3 

27 

40 

22 

61 

70 
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Table 6. Obtained Results for the 33-bus RTS for Three DG Units 
Location Size (P [MW]) Performance evaluation indices 


CASE 

# 

DGi 

dg 2 

DG, 

DG, 

dg 2 

DG, 

PF 

DGpl 

^max 

P Losses 

VSI 

VP 

ALR 

RLR 

SLI 

VPI 

VSIM 

NBVV 

CASE 

1.1 

- 

- 

- 

- 

- 

- 

- 

- 

3.406 

0.211 

1 

0.134 

- 

- 

- 

- 

- 

21 

CASE 

1.22 

14 

17 

32 

0.310 

0.310 

0.311 

1 

25 

3.843 

0.118 

1 

0.053 

43.865 

45.172 

12.814 

60.272 

0 

5 

CASE 

1.23 

15 

18 

32 

0.618 

0.617 

0.619 

1 

50 

4.187 

0.107 

1 

0.015 

49.263 

47.059 

22.920 

88.422 

0 

0 

CASE 

1.24 

15 

17 

32 

0.928 

0.929 

0.929 

1 

75 

4.463 

0.155 

0.865 

0.010 

26.562 

17.411 

31.004 

92.732 

13.511 

0 

CASE 

1.25 

10 

16 

32 

1.238 

1.238 

1.238 

1 

100 

4.684 

0.189 

0.828 

0.016 

10.421 

0.843 

37.497 

87.801 

17.243 

0 

CASE 

1.26 

14 

17 

31 

1.003 

0.746 

2.926 

1 

- 

4.977 

0.296 

0.823 

0.025 

40.088 

57.031 

46.098 

80.953 

17.698 

0 

CASE 

1.27 

10 

14 

32 

0.296 

0.296 

0.296 

0.814 

25 

3.916 

0.083 

1 

0.035 

60.522 

62.071 

14.965 

73.469 

0 

0 

CASE 

1.28 

14 

17 

32 

0.627 

0.630 

0.630 

0.865 

50 

4.473 

0.059 

0.867 

0.010 

71.836 

68.994 

31.305 

92.621 

13.344 

0 

CASE 

1.29 

18 

31 

32 

0.912 

0.910 

0.903 

0.835 

75 

4.861 

0.077 

0.828 

0.014 

63.663 

53.814 

42.700 

89.359 

17.187 

0 

CASE 

1.30 

15 

30 

31 

1.232 

1.151 

1.192 

0.905 

100 

5.052 

0.097 

0.823 

0.029 

53.882 

47.934 

48.291 

78.095 

17.710 

0 

CASE 

1 

14 

17 

30 

0.906 

0.185 

2.212 

0.810 

- 

5.063 

0.091 

0.823 

0.032 

56.902 

51.934 

48.635 

75.889 

17.724 

0 



Figure 5. System maximum loading curves for CASE 1.1, CASE 1.5 and CASE 1.10 
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Figure 6. Voltage profiles for CASE 1.1, CASE 1.5 and CASE 1.10 
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Figure 7. Evolution of the objective function for CASE 1.5 and CASE 1.10 
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Figure 8. System maximum loading curves for CASE 1.1, CASE 1.15 and CASE 1.20 



Figure 9. Voltage profiles for CASE 1.1, CASE 1.15 and CASE 1.20. 
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Figure 10. Evolution of the objective function for CASE 1.15 and CASE 1.20 



Figure 11. System maximum loading curves for CASE 1.1, CASE 1.25 and CASE 1.30 
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Figure 12. Voltage profiles for CASE 1.1, CASE 1.25 and CASE 1.30 
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IttiMioni 

Figure 13. Evolution of the objective function for CASE 1.25 and CASE 1.30 


Results for the 69-bus radial test system 

Results for one DG unit incorporation 

The obtained results for the incorporation of one DG unit for the 69-bus RTS are given in Table7. 

The following comments can be made on this table: 

a. For the base case without DG (i.e. CASE 2.1) the maximum loadability is 3.21. In this case NBVV=9. 

b. For CASE 2.2 NBVV is equal to 3. Even though the NBVV has been decreased from 9 in the base case 
to 3 in this case, the network is still not fully optimized due to incorporation of only one DG and the 
limitation imposed on the penetration level (25% only). 

c. DG units are located mainly in busses 61, 62, 63 and 64. 

d. For CASE 2.3, using only one DG unit but with a penetration level of 50% the network is optimized and 
the NBVV=0. The A max in this case has been increased from 3.210 to 3.966 compared with the base 
case i.e. SLI=23.55%. Likewise, for CASE 2.4 and CASE 2.5 with DG PL =75% and DG PL =100%, the 
SLI is equal to 34.22% and 44.86%, respectively. 

e. For CASE 2.7, even though only one DG unit is incorporated using a penetration lever of 25% the 
network is optimized (NBVV=0) compared with CASE 2.2 because the PF has been computed as 0.89 
instead of 1 for CASE 2.2. 

f. The PV curves of CASE 2.1, CASE 2.5 and CASE 2.10 are displayed in Figure 14. 

g. In all cases, voltage profile is better than the base case as shown by the VPI. The best improvement is 
obtained for CASE 2.4 where VPI=88.14%. The voltage profile of CASE 2.1, CASE 2.5 and CASE 
2.10 are displayed in Figure 15. 

h. The active power losses is reduced in all cases. 

i. The evolution of the objective function for CASE 2.5 and CASE 2.10 is given in Figure 16. 

Results for two DG units incorporation 

The obtained results for the incorporation of two DG units for the 69-bus RTS are given in Table 8. 

The following comments can be made on this table: 

a. When incorporating two DG units only one case is not optimized where NBVV=3 that is CASE 2.12. 

b. The first DG unit is located at buses 38, 61, 62, and 63 while the location of the second DG unit is 
mainly on buses 61, 62, 63, 64 and 65. 

c. The most important SLI is for CASE 2.21 where it is equal to 53.25%. An illustration of the PV curves 
of CASE 2.1, CASE 2.15 and CASE 2.20 is given in Figure 17. 

d. The voltage profile is improved in all cases compared with the base case. An illustration of VP for 
CASE 2.1, CASE 2.15 and CASE 2.20 is given in Figure 18. 

e. Active power losses are reduced in all cases. 

f. When the PF is not fixed and it is considered as a design variable, it varies approximatively between 
0.80 and 0.90. 

g. The evolution of the objective function for CASE 2.15 and CASE 2.20 is sketched in Figure 19. 

Results for three DG units incorporation 

The obtained results for the incorporation of three DG units for the 69-bus RTS are shown in 

Table 9. The following comments can be made on this table: 

a. When incorporating three DG units only one case is not optimized where NBVV=4 that is CASE 2.22. 
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b. The best SLI obtained is for CASE 2.31 where it is equal to 53.11%. An illustration of the PV curves of 
CASE 2.1, CASE 2.21 and CASE 2.25 is given in Figure 20. 

c. The voltage profile is improved in all cases compared with the base case. An illustration of VP for 
CASE 2.1, CASE 2.21 and CASE 2.25 is given in Figure 21. 

d. Active power losses are reduced in all cases compared with the base case. 

e. When the PF is not fixed and it is considered as a design variable, it varies approximately between 0.81 
and 0.98. 

The evolution of the objective function for CASE 2.21 and CASE 2.25 is given in Figure 22. 


Table 7. Obtained Results for the 69-bus RTS for One DG Unit 


CASE Location Size (P [MW]) Performance evaluation indices 


# 

DGi 

dg 2 

DG, 

DG, 

dg 2 

DG, 

PF 

DGpl 

^max 

P Losses 

VSI 

VP 

ALR 

RLR 

SLI 

VPI 

VSIM 

NBVV 

CASE 

2.1 

- 

- 

- 

- 

- 

- 

- 

- 

3.210 

0.225 

1 

0.099 

- 

- 

- 

- 

- 

9 

CASE 

2.2 

64 

- 

- 

0.951 

- 

- 

1 

25 

3.604 

0.118 

1 

0.046 

47.736 

44.796 

12.275 

53.952 

0 

3 

CASE 

2.3 

64 

- 

- 

1.900 

- 

- 

1 

50 

3.966 

0.099 

1 

0.019 

55.938 

52.470 

23.546 

80.803 

0 

0 

CASE 

2.4 

62 

- 

- 

2.851 

- 

- 

1 

75 

4.309 

0.118 

0.940 

0.012 

47.456 

47.201 

34.223 

88.174 

5.969 

0 

CASE 

2.5 

61 

- 

- 

3.801 

- 

- 

1 

100 

4.650 

0.198 

0.843 

0.020 

11.904 

16.461 

44.855 

80.339 

15.692 

0 

CASE 

2.6 

61 

- 

- 

4.011 

- 

- 

1 

200 

4.725 

0.223 

0.823 

0.023 

0.919 

6.733 

47.182 

76.688 

17.689 

0 

CASE 

2.7 

64 

- 

- 

1.029 

- 

- 

0.885 

25 

3.771 

0.071 

1 

0.033 

68.505 

64.326 

17.477 

66.536 

0 

0 

CASE 

2.8 

64 

- 

- 

2.011 

- 

- 

0.863 

50 

4.307 

0.046 

0.938 

0.012 

79.370 

74.578 

34.178 

88.095 

6.156 

0 

CASE 

2.9 

63 

- 

- 

2.999 

- 

- 

0.858 

75 

4.817 

0.085 

0.827 

0.020 

62.283 

61.809 

50.057 

79.611 

17.266 

0 

CASE 

2.10 

61 

- 

- 

3.008 

- 

- 

0.800 

100 

4.920 

0.093 

0.823 

0.022 

58.603 

60.160 

53.251 

77.565 

17.712 

0 

CASE 

7 1 1 

61 

- 

- 

3.008 

- 

- 

0.800 

200 

4.920 

0.093 

0.823 

0.022 

58.603 

60.160 

53.251 

77.565 

17.712 

0 


Table 8. Obtained Results for the 69-bus RTS for Two DG Units 


Location Size (P [MW]) Performance evaluation indices 


# 

DG, 

dg 2 

DG, 

DG, 

dg 2 

DG, 

PF 

DGpl 

^max 

PLosses 

VSI 

VP 

ALR 

RLR 

SLI 

VPI 

VSIM 

NBVV 

CASE 

2.1 

- 

- 

- 

- 

- 

- 

- 

- 

3.2 

10 

0.225 

1 

0.0 

99 

- 

- 

- 

- 

- 

9 

CASE 

63 

65 


0.47 

0.475 


1 

25 

3.6 

0.117 

1 

0.0 

48.1 

45.2 

12.1 

53.8 

0 

i 

2.12 


5 


1 

00 

46 

21 

20 

38 

74 

j 

CASE 

62 

64 


0.95 

0.949 


1 

50 

3.9 

0.087 

1 

0.0 

61.3 

58.5 

23.7 

80.7 

0 

0 

2.13 


0 


72 

19 

93 

61 

28 

85 

CASE 

2.14 

61 

64 

- 

1.42 

4 

1.424 

- 

1 

75 

4.3 

29 

0.121 

0.925 

0.0 

12 

46.1 

58 

45.7 

10 

34.8 

62 

87.7 

59 

7.540 

0 

CASE 

2.15 

61 

63 

- 

1.90 

0 

1.901 

- 

1 

100 

4.6 

57 

0.201 

0.836 

0.0 

20 

10.4 

87 

14.8 

47 

45.0 

83 

79.5 

92 

16.377 

0 

CASE 

2.16 

61 

65 

- 

3.96 

7 

0.044 

- 

1 

200 

4.7 

26 

0.223 

0.823 

0.0 

23 

1.00 

8 

6.82 

7 

47.2 

28 

76.6 

13 

17.686 

0 

CASE 

2.17 

63 

64 

- 

0.52 

3 

0.523 

- 

0.8 

98 

25 

3.7 

65 

0.068 

1 

0.0 

33 

69.5 

76 

65.5 

70 

17.2 

94 

66.3 

08 

0 

0 

CASE 

2.18 

63 

64 

- 

0.95 

2 

0.952 

- 

0.8 

17 

50 

4.3 

07 

0.031 

0.962 

0.0 

11 

86.1 

19 

82.0 

68 

34.1 

78 
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43 
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0 

CASE 

2.19 

61 

64 

- 

1.40 

9 

1.410 

- 

0.8 

08 

75 
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44 
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20 

64.5 

05 
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65 
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78 
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01 
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0 
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62 

- 
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1 
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- 

0.8 

06 
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14 
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22 

59.6 

12 
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88 
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69 
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30 

17.727 

0 

CASE 

2.21 

38 

61 

- 
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2 

3.008 

- 
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00 
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4.9 

20 

0.093 
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22 
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30 
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30 
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51 
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76 
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Table 9. Obtained Results for the 69-bus RTS for Three DG Units 
CASE Location Size (P [MW]) Performance evaluation indices 


# 

DG, 

dg 2 

dg 3 

DG, 

dg 2 

dg 3 

PF 

DGpl 

^•max 

P Losses 

VSI 

VP 

ALR 

RLR 

SLI 

VPI 

VSIM 

NBVV 

CASE 

2.1 

- 

- 

- 

- 

- 

- 

- 

- 

3.210 

0.225 

1 

0.099 

- 

- 

- 

- 

- 

9 

CASE 

2.22 

61 

62 

65 

0.317 

0.317 

0.317 

1 

25 

3.597 

0.115 

1 

0.046 

48.915 

46.094 

12.047 

53.453 

0 

4 

CASE 

2.23 

61 

63 

64 

0.633 

0.631 

0.633 

1 

50 

3.969 

0.085 

1 

0.019 

62.313 

59.584 

23.637 

80.675 

0 

0 

CASb 

2.24 

17 

62 

64 

0.806 

0.948 

0.949 

1 

75 

4.011 

0.079 

0.990 

0.002 

64.887 

62.013 

24.960 

97.994 

0.994 

0 

CASE 

2.25 

61 

62 

64 

1.267 

1.267 

1.266 

1 

100 

4.674 

0.205 

0.826 

0.021 

9.094 

13.241 

45.585 

78.702 

17.429 

0 

CASE 

2.26 

17 

46 

61 

0.008 

0.023 

4.010 

1 

200 

4.725 

0.223 

0.823 

0.023 

1.064 

6.860 

47.182 

76.850 

17.693 

0 

CASb 

2.27 

59 

63 

65 

0.340 

0.338 

0.339 

0.876 

25 

3.713 

0.072 

1 

0.035 

68.066 

64.199 

15.651 

64.794 

0 

0 

CASE 

2.28 

63 

64 

65 

0.757 

0.757 

0.757 

0.975 

50 

4.249 

0.065 

0.937 

0.012 

70.901 

67.222 

32.352 

87.522 

6.272 

0 

CASE 

2.29 

9 

61 

64 

0.802 

0.979 

0.979 

0.841 

75 

4.357 

0.021 

0.943 

0.007 

90.507 

87.183 

35.729 

93.240 

5.695 

0 

CASb 

2.30 

62 

64 

67 

1.449 

1.451 

1.444 

0.934 

100 

4.725 

0.081 

0.824 

0.017 

64.029 

61.687 

47.182 

83.358 

17.574 

0 

CASb 

7 

1 

32 

61 

2.266 

0.885 

3.029 

0.809 

200 

4.915 

0.097 

0.823 

0.022 

56.749 

58.053 

53.114 

77.415 

17.698 

0 
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Figure 14. System maximum loading curves for CASE 2.1, CASE 2.5 and CASE 2.10 
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Figure 15. Voltage profiles for CASE 2.1, CASE 2.5 and CASE 2.10 
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Figure 16. Evolution of the objective function for CASE 2.5 and CASE 2.10 
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Figure 17. System maximum loading curves for CASE 2.1, CASE 2.15 and CASE 2.20 
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Figure 18. Voltage profiles for CASE 2.1, CASE 2.15 and CASE 2.20 
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Figure 19. Evolution of the objective function for CASE 2.15 and CASE 2.20 



Figure 20. System maximum loading curves for CASE 2.1, CASE 2.21 and CASE 2.25 
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Figure 21. Voltage profiles for CASE 2.1, CASE 2.21 and CASE 2.25 
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Iterations 

Figure 22. Evolution of the objective function for CASE 2.21 and CASE 2.25 


5. CONCLUSION 

This paper has investigated the optimal simultaneous incorporation of multiple distributed 
generation units for the maximization of the system loadability. The maximum loadability of the system is 
calculated using an accurate and fast approach compared with classical approach. The time gained by the 
proposed approach can be huge in an optimization procedure as for OIDG. 

The OIDG in this paper has been solved using the vortex search algorithm. Many cases have been 
investigated using one, two and three DG units. Moreover, different penetration level of DG have been 
considered. In addition, in some cases the power factor is considered as unity and some other cases it is 
considered as a design variable in order to maximize the system loadability. From the obtained results, it can 
be concluded that; system loadability depends upon the level of DG penetration level. System loadability 
improves when the power factor of the generator changes i.e. when the generator supplies reactive power 
also. When the number of DG units is increased, the maximum loadability also improves. When the number 
of DG units is increased, the power system losses are reduced. Total number of DG units and their power 
factor changes the system loadability as well as power factor. Future research should concentrate on the 
investigation of a multi-objective OIDG problem. 
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